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We consider a magnetic flux pointing in the 2 direction of an axially symmetric space-time (Melvin 
Universe) in a Born-Infeld-type extension of General Relativity (GR) formulated in the Palatini 
approach. Large magnetic fields could have been produced in the early Universe, and given rise 
to interesting phenomenology regarding wormholes and black hole remnants. We find a formal 
analytic solution to this problem that recovers the GR result in the appropriate limits. Our results 
set the basis for further extensions that could allow the embedding of pairs of black hole remnants 
in geometries with intense magnetic fields. 
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I. INTRODUCTION 

The Melvin Universe is a regular, non-black hole elec¬ 
trovacuum solution of the Einstein-Maxwell equations 
describing a bundle of magnetic flux lines in static equi¬ 
librium, held together by their own gravitational in¬ 
teraction p3]. It represents an explicit realization of 
Wheeler’s geons [2], namely, self-consistent sourceless 
gravitational-electromagnetic entities Q. This solution 
has been generalized to rotating and time-dependent con¬ 
figurations H, dilatonic [H] and axion fields @, higher 
dimensions p}, and nonlinear electrodynamics Q- Its 
applications include supergravity jf)§ and D-Branes 0- 

The magnetic background of the Melvin solution can 
be employed to investigate the creation of black hole 
pairs, since the negative energy of the magnetic field com¬ 
pensates the energy of the pair and thus energy conser¬ 
vation is fulfilled [llj . If the magnetic field is strong 
enough, a Wheeler wormhole solution, namely, a pair 
of extreme Reissner-Nordstrom black holes identified at 
their throats and with opposite charges, can be formed 
0] ■ This construction faces us with the problematic 
issue of topology change processes, in which the quan¬ 
tum foam scenarios are based. In the case of the Melvin 
space-time, such a process is governed by tunelling ef¬ 
fects through instantons, which can be modeled in terms 
of Ernst’s metric fl3| . 

Large-scale coherent magnetic fields with strength of 
order 10 -6 G are observed in galaxies and galaxy clus¬ 
ters 0 . According to the primordial hypothesis, these 
magnetic fields were created in the early Universe and 
they were later amplified by a dynamo mechanism (see 
e.g. fl5l] for a review). Production of microscopic black 
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holes/wormholes in the early Universe by large fluctu¬ 
ations was first discussed by Hawking 0 ], with the re¬ 
sult that large numbers of such objects with Planck-order 
mass and a few units of charge could have been produced. 
Though the interest on such an idea faded away after 
the discovery of Hawking radiation, which would imply 
that such objects should have evaporated by the current 
epoch, recent research in extensions of GR using different 
approaches has put again regular black holes, wormholes, 
and black hole remnants under an intense discussion (see 
0 for a review). Indeed, a number of scenarios have 
been proposed 0 and some of them have the necessary 
ingredients for the production of microwormholes. Astro- 
physical observations put some constraints on the pres¬ 
ence of wormholes in our Universe, but their existence 
cannot be ruled out 0 ]. For instance, some traversable 
wormholes are viable candidates to explain the supermas- 
sive objects at the centers of galaxies, while other kinds 
of wormholes can be excluded 0. 

The main goal of this paper is to work out a Melvin- 
type space-time in an extension of GR which has gained 
interest in the last few years, namely, Born-Infeld grav¬ 
ity. Recently, it was shown [21], 0] that this theory con¬ 
tains static, spherically symmetric wormholes supported 
by the electromagnetic field with geonic properties. A 
very important question is to discuss plausible mecha¬ 
nisms for the generation of such wormholes. An impor¬ 
tant step in this sense was done in 0], where it was 
shown that dynamical generation of wormholes through 
charged fluxes of radiation is possible. In this work we 
progress in a different, though related, direction hoping 
to clarify if such wormholes could be generated (or em¬ 
bedded) in highly magnetized scenarios. The first step 
requires to determine if Melvin-type solutions in analyt¬ 
ical form can be found in this framework. We find that 
the answer is partially positive, as two analytical approx¬ 
imations to the exact solution can be found explicitly in 
the two regimes of interest. The fact of having some an¬ 
alytical control on the solutions suggests that an answer 
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to the problem posed here might be accessible, though 
further research is still necessary. 

The paper is organized as follows: in Sec. m we briefly 
review Born-Infeld gravity in the Palatini approach. In 
Sec. EH we work out the axially symmetric Melvin-type 
space-time in this scenario, and obtain analytical solu¬ 
tions in the two regions of interest, namely, close to the 
axis and far from it. Conclusions and future research are 
discussed in Sec. m 


II. ACTION AND MAIN EQUATIONS 


In the last few years an extension of GR, that has at¬ 
tracted much attention in astrophysical and cosmological 
scenarios [24j ]. has been proposed following the analogy 
with the Born-Infeld theory of non-linear electrodynam¬ 
ics [ 25 } . Initially introduced in the metric formalism [26| , 
the interest in this theory arose once its Palatini version 
was considered as it avoids higher-order derivatives 
and ghosts. The action of this Born-Infeld theory of grav¬ 
ity (BI for short) coupled to an electromagnetic field can 
be written as 

s BI = ^- [fx^-W^g] 

k e j 

- . (1) 

where n 2 is Newton’s constant in suitable units (in GR, 
n 2 = 87tGn), g and q are the determinant of the space- 
time metric g^ and of the metric q^ = g^u + eR^, with 
e a small constant with dimensions of length squared, 
and = R a fj.au is the (symmetric) Ricci tensor 

constructed with the affine connection T = r^„, which 
is a priori independent of the metric g^u (metric-affine 
or Palatini approach). In the matter (Maxwell) sector, 
Ffiu = dfjAfj — d„Afj is the field strength tensor of the 
vector potential The meaning of A in Eq. (JT]) follows 
from a series expansion in powers of e of the gravitational 
sector, Sq, as 


Sg 


d A Xy/~g 
2 K 2 


R-2A-| 



+ RfjuR^ 



to support both the Palatini approach and the Born- 
Infeld gravity action j2f|, with potential consequences 
for our understanding of the microscopic description of 
space-time and gravitational phenomena [30 ]. 

To implement the Palatini formalism, we perform in¬ 
dependent variations of the action © with respect to the 
metric and the connection, which yields (setting torsion 
to zero for simplicity [El) 

qi ju _ Xgf ju = _ K 2 eT ^ j ( 2 ) 

V^(v^r)=0. (3) 

From Eq. ©, it follows that the independent connection 
is given by the Christoffel symbols of the metric q^. 
Using 0, one finds that q^u is related to the space-time 
metric g M[/ (for which one has V^( \/—gg^ u ) ^ 0) as 


VW\ 

V\i\ 



( 4 ) 


with the definition 


E = XI- en 2 f , 


( 5 ) 


where T = T^ a g a u and a hat denotes a matrix. This 
clearly shows that the relation between q^ and is 
algebraic and only governed by the matter fields. 

From the definition q = g + eR and the relations 0, 
a bit of algebra allows us to write the metric field equa¬ 
tions © in terms of q^ as 


V(?) 



(£ g <5/ + V) , 


( 6 ) 


where Lq = (|S| 1//2 — A)/(e/« 2 ) is the gravity Lagrangian. 
This represents a set of second-order differential Einstein- 
like field equations with all the right-hand-side only de¬ 
pending on the matter sources. Since q„„ is algebraically 
related to g M1/ , the field equations for it are second-order 
as well. In vacuum, the field equations m boil down to 


Rfusio.) 


(A-l) 

Ae 


Qfjv Rfjuifj) 


(A — 1) 


9^ ■ 


This theory recovers GR with a cosmological constant 
term A = ^ at zeroth order, while higher-order correc¬ 

tions in the curvature invariants are suppressed by powers 
of e. 

In the Palatini formulation of modified gravity, the 
usual troubles with higher-order field equations and 
ghost-like instabilities are avoided for large families of 
models. This is in sharp contrast with the situation 
in most approaches to modified gravity, where the con¬ 
nection is taken to be a priori given by the Christoffel 
symbols of the metric. While in the special case of the 
Einstein-Hilbert action both approaches give the same 
Einstein equations, this is not so for most extensions of 
GR [28] . We point out that some experimental results re¬ 
garding systems with defects in solid state physics seem 


This clearly shows that the dynamics of this theory in 
vacuum is that of GR with a cosmological constant term 
with value A = iAzll j n addition, it implies that the 
theory is free of extra propagating degrees of freedom 
and ghost-like instabilities. 

III. AXIAL MAGNETIC SPACE-TIME 
A. Basic equations 

The electromagnetic (Maxwell) field in action © sat¬ 
isfies the equations 

= 0 dfj (V=gF>”) = 0 . 


( 7 ) 











3 


For a radial, static, and spherically symmetric field, the 
only non-vanishing component is F tr = E(r). In this 
case, one finds that if e = — 21 %, with l e representing 
a length scale, then the point-like singularity of GR 
is generic.ally replaced by a wormhole structure, whose 
properties have been studied in detail in a number of pa¬ 
pers [U [22j . We will thus adopt from now on the choice 
e = —21%. Our interest here is to find the analog solution 
for the case of cylindrically symmetric and magnetically 
charged configurations in coordinates (t,z,p,<j>). The 
physical interpretation is that of a space-time sourced by 
a beam of magnetic fields parallel to the z-axis. In order 
to solve the system of field equations (0 for BI grav¬ 
ity o with an axially symmetric electromagnetic field 
described by Maxwell equations 0, we introduce a line 
element for the metric g^ as follows 


gravitational length scale l e to the magnetic length scale 
If 3 controls the departures of our solutions from those of 
GR. 

Since the field equations 0 will be solved in terms 
of the metric q pl , associated to the affine structure, we 
define an axially symmetric line element for it as 

ds 2 = f(p)(—dt 2 -I- dz 2 ) + g(p)dp 2 + h{p)p 2 d <\> 2 , (14) 

which is formally identical to that of the metric g [see 
Eq. ©]. From the relations (1121) . we immediately see 
that the functions defining the line elements of the metric 
(El) and g pv ® are related by 

/ = cr+/, h = a-h, g = <r-g. (15) 

From the first of the above relations, it is easy to see that 


ds 2 = f(p)(—dt 2 + dz 2 ) + g(p)dp 2 + h(p)p 2 dtp 2 . (8) 

With this line element, from the field equations 0 we 
find that the only non-vanishing component of the field 
strength tensor is F pv , which satisfies 

8 ,(^-) = 0 = ( 9 ) 

where /3 is an integration constant that determines the 
intensity of the magnetic field. With this result, one can 
verify that the field invariant is 7),,, F IJV = 2 F ptp F PV = 
2 (F pip ) 2 ghp 2 . Accordingly, to get rid of the numerical 
factor, we define X = — \F pu F fJjl ' = —/3 2 /f 2 . We note 
that X depends neither on g nor on h , which simplifies 
the analysis. With these results, we find that the energy- 
momentum tensor of the electromagnetic field 

V = ~ (io) 


/ + \Jf 2 -4f? 
2 


(16) 


which indicates that / is bounded by / > 2 f c . For this 
value of /, we get / = / c , which corresponds to er_ = 0. 


B. Computation of the metric components and the 
line element 


With the elements above, we note that the BI La- 
grangian can be written as 


V det E - 1 _ /? 2 / 2 
— 2 k 2 1 % 87r/ 4 


(17) 


It is a simple calculation to show that the field equa¬ 
tions 0 for these axially symmetric solutions take the 
explicit form 


takes the simple form 

T/ = |- d iag(l,l,-l,-l) . (11) 

The matrix E of the theory, E = 7 + 2 1%k 2 T, then be¬ 
comes 


(7—i 6 \ 

6 a + i)' 


( 12 ) 


where we have defined 


<t± = 1±£ (13) 

To simplify the notation, here we have denoted / 2 = 
k 2 /3 2 l 2 , with k 2 = k 2 /4:Tt , and 0 and 7 are the zero and 
identity 2x2 matrices, respectively. Note that k 2 /3 2 can 
be interpreted as an inverse squared length associated to 
the magnetic field, which we denote l 2 p = l/k 2 /3 2 . This 
leads to f c = l e /lp and indicates that the ratio of the 




d 2 0 2 ( ^ 6 \ 

2 f 2 Vo -±i ) • 


(18) 


Using the xAct package for Mathematica HH, the com¬ 
ponents of 7? /i !/ (g) corresponding to the line element (fill) 
can be easily obtained, and are given by 


7?o° = i?i 1 =- 


R2 2 = 


4/2 




2/1 


pp 


(19) 


£Z 

Ah 2 p-j- + h (pfphp + 2 h{f p - 2 pfpp^j 


4 f 2 h 2 p 

+ f (ph 2 p ~ 2h(2h p + phpp )) 

R 3 13 = . -„V„ \f{ph 2 p - 2h(2h p + phpp )) 


4 f 2 h 2 p 
h(pfphp + 2 hf^j 


( 20 ) 


( 21 ) 


where we have imposed the gauge freedom g = /, and 
have used a subindex to denote derivative with respect 
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to p. Given the block-diagonal form of the field equa¬ 
tions (USD, it follows that i ? 2 2 — -R 3 3 = 0 , which implies 
that 


0 = 


1 


4 f 2 h 2 p 


f 2 

4 ~ h2 P~y + ~ h ( 2 Pfph P + 4h{/ p - p/pp}) 


( 22 ) 

With elementary algebraic manipulations, this equation 
becomes 


(hp 2 2/ p \ _ 2/pp 

U p /) h ' 

which can be readily integrated to obtain 



(23) 


(24) 


where a is an integration constant (with dimensions), 
whose value we keep undefined for the moment. 

Using now K3l) in the expression for R 0 °, we get the 
following equation for /: 


f 3 /p « 2 /? 2 / 

Jpp 4 f 2 / ' 


(25) 


We will next proceed to solve this equation to obtain 
f(p). With the solution, the geometry g p„ can be ob¬ 
tained using the relations fl]). 


C. Finding a solution 


Using the expression (flfil) in (l25l) . the equation to solve 
takes the form 


? 3 ft _ 

Ipp 4 f ~ - r~ - 

4/ / + V/ 2 ^ 4 /c 


(26) 


This equation can be written in a more convenient (di¬ 
mensionless) form by defining / = 2 f c (j){x) and p 2 = 
which leads to 


_ = 0 
4 0 0 + ' 


(27) 


We now define a new function 12 = 0 2 such that dfl/d(f> = 
2 (f> xx , which turns (EH) into 


where C is an integration constant and 2 U a hypergeo¬ 
metric function. Given that 12 = 0 2 must be a positive 
function by construction, an expansion about 0 « 1 in¬ 
dicates that C cannot be arbitrary. In fact, we find that 


lim 12 ss C + — 
<£—>i 3 


V^(f) 

r(D 


+ (<£-!) 



4V^r(§) \ 

r (!) J 


(30) 

+ 0 [(^-l) 2 ] , 


which requires that the constant factor when 0 = 1 be 
positive or zero. Now, given that 0 must be greater 
or equal to 1 to guarantee the reality of the differential 
equation (l27l) , on physical grounds we must demand that 
0 = 1 be a minimum. This forces us to set this constant 

to zero, which implies C = — | + « 2.16274. We 

note that in GR the constant C is usually normalized as 
Cgr = 2. 

The expression for 0(x) can be obtained numerically 
by integrating <p x = fl 1 / 2 with the value of C chosen 
above. It is useful, however, to find approximate analyti¬ 
cal solutions in the two regimes of interest, namely when 
0 —>■ 1 and when 0 1. The limit 0 1 is particu¬ 

larly interesting because the equation for 0 (x) converges 
to the equation that one finds in GR. In fact, in this limit 
m turns into 


(31) 

and proceeding as above, we get flcR — Ccf. > 3//2 — 20 
(where C 0 3 / 2 represents the solution of the homoge¬ 
neous equation). With this solution one can readily find 
0(x) = 4(1 + (Cx) 2 /32) 2 jC 2 , which is the right GR 
Melvin solution Jj. The choice Cgr = 2 arises natu¬ 
rally and makes the metric Minkowskian as x —> 0, where 
0(x) —> 1. In our case, however, C ^ 2 is motivated by 
the behavior of the differential equation near <f> —> 1. In 
this other limit, we find that 

<t>l « 2(0 - 1) , (32) 

which leads to 0(x) ~ 1 + x 2 /2. With this result, one 
can verify that the line element as x —> 0 takes the form 


ds 2 

T 


(l+a 


— dt 2 + dz 2 + pi 


dx 2 


2 fcPl 


xdip 2 , (33) 


12^ 



20 

0 + 7^1 


(28) 


This equation admits an exact solution for 12(0) of the 
form 


f2 = G0l + ^ (0-70^1) 

_ 80 (i i .5 n 

3 2 Fl U’ 2 ’4’ 0 ^ > 


( 29 ) 


where p$ = 2 l e lp and a is the (unspecified) integration 
constant that arose in (UHl) . Under a rescaling of the x 
coordinate of the form dx 2 lx = dy , this line element 
becomes 


ds 2 

T 


- dt 2 + dz 2 + ^ dy 2 


+ y 2 dp 2 , (34) 
°JcP o 


where second-order corrections in y 2 have been neglected. 
This last expression shows that the Minkowski space-time 























5 


^Num 

^App 



Figure 1. Ratio of the function f(x) computed numerically 
by its analytic approximation used in (13511 . 


can be recovered near the axis by just defining r = poy/2 
and a = 2/ c po- An additional global rescaling of units 
could be used to absorb the constant factor f c . 

In the (j) 1 limit, we find that / « 2 f c <f>. Taking the 

solution for <f> obtained above in this limit, (f>(x ) = 4(1 + 
(Cx) 1 2 3 4 /32) 2 / C 2 , the line element can be approximated as 



where we have taken a = 2/ c /9q as above. In Figjl] we 
show the ratio between the function /( x) computed nu¬ 
merically and its approximation for 1. The ap¬ 

proximation is very good for values of x > 20. It is 
worth noting that the line element (1551) can be made 
to agree with the GR solution (up to a constant con¬ 
formal factor) by just introducing a constant rescaling 
of ( t,z,p ) —>■ ( Xt,Xz,Xp ) with A 2 = 64 f c /C 2 , and by 
suitably choosing the integration constant a in the cor¬ 
responding definition of hp 2 . 


IV. SUMMARY AND PERSPECTIVES 

In this work we have investigated the possibility of find¬ 
ing analytical solutions for axially symmetric magnetic 
fields in the Born-Infeld theory of gravity. This type of 
non-asymptotically flat solutions are known as Melvin 
universes and are of great interest to give plausibility to 
the generation of pairs of entangled black holes by intense 
magnetic fields. In the context of GR these solutions are 
well known and there exist advanced solution-generating 
methods that allow to embed electrically charged solu¬ 
tions within these magnetized scenarios in a very elegant 
and robust way, even in the case of stationary space- 
times. For the Born-Infeld theory of gravity there is no 
guarantee that such methods can be implemented or even 
exist at all. For this reason, as a first step in this direc¬ 
tion, we have investigated the very existence of Melvin- 
type solutions. By a suitable choice of variables (using 
an auxiliary line element), we have been able to write 
the field equations in a simple dimensionless form and 
obtain analytical solutions in the two regimes of interest, 
namely, near the symmetry axis and far from it. The 
complete solution can be easily worked out numerically 
and we have shown that our analytical approximations 
fit well with the exact computation. As a result, near 
the axis one recovers a Minkowskian geometry and far 
from it a standard Melvin Universe. This suggests that 
the embedding of electrically charged solutions, includ¬ 
ing pairs of black hole remnants and wormholes, might 
be possible at least in some approximate form, which mo¬ 
tivates further research in this direction. To investigate 
these questions in detail one should consider the analog 
of Ernst metric [l3| in our scenario. Progress in this sense 
is currently underway. 
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